Casimir energy for a scalar field with a frequency dependent boundary condition 



H. Falomir, K. Rebora 

IFLP - Dept. de Ftsica, Fac, de Ciencias Exactas, Universidad Nacional de La Plata, C.C. 67, (1900) La Plata, Argentina 

and M. Loewe 

Facultad de Ftsica, Pontificia Universidad Catolica de Chile, Casilla 306, Santiago 22, Chile 

(February 1, 2008) 



O 

o 
o 

(N 

> 

O 

O 



> 

in 

oo 

o 

O 
O 



Oh 



X 



We consider the vacuum energy for a scalar field subject 
to a frequency dependent boundary condition. The effect of 
a frequency cut-off is described in terms of an incomplete (- 
function. The use of the Debye asymptotic expansion for 
Bessel functions allows to determine the dominant (volume, 
area, . . . ) terms in the Casimir energy. The possible interest 
of this kind of models for dielectric media (and its application 
to sonoluminescence) is also discussed. 
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I. INTRODUCTION 

The Casimir effect arises as a distortion of the 
vacuum energy of quantized fields due to the presence of 
boundaries (or nontrivial topologies) in the quantization 
domain. This effect, which has a quantum nature associ- 
ated with the zero-point oscillations in the vacuum state, 
is significant in diverse areas of physics, from statistical 
physics to elementary particle physics and cosmology. 

In particular, in the last years there has been a great 
interest in the Casimir energy of electromagnetic fields 
in the presence of dielectric media, due to Schwinger's 
suggestion that it could play a role in the explanation 
of the phenomenon of sonoluminescence Q . 

The results obtained on this subject by different groups 
through several calculation techniques (as Green's func- 
tions methods, van der Waals forces, ^-function methods 
and asymptotic developments for the density of states - 
see references |p|-p^ among others) are rather controver- 
sial, and some basic issues remain to be clarified. In this 
respect, it is our aim to contribute to the understand- 
ing of the problem by studying a simplified model, which 
incorporates a frequency cut-off in the boundary condi- 
tions at the separation between media, to emulate the 
behavior of real dielectrics. 

It should be mentioned that authors introducing a cut- 
of in the wave number to describe the behavior of real 
dielectrics agree with Schwinger's explanation of the phe- 
nomenon (See, for example, |lT|-p^). On the other hand 
in Ref. p5| , where the permittivity of a non-magnetic 
medium is modeled by means of a (angular momentum 
dependent) one-absorption-frequency Sellmeir relation, 
and a (^-function technique is employed to sum-up the 



contributions of the proper oscillation modes of the elec- 
tromagnetic field, it is concluded that the vacuum energy 
of a bubble embedded in this material has the wrong be- 
havior with the radius (besides its absolute value being 
far too small) to be relevant to sonoluminescence. 

Candelas has been the first to remark the impor- 
tance of dispersion in connection with the Casimir ef- 
fect. Later, Brevik and collaborators p7|-[2^ have con- 
sidered several models with spherical and cylindrical ge- 
ometries, showing that dispersive effects may reverse the 
direction of the Casimir force acting on the boundary 
between media. In Ref. it is studied the case 

of spherically symmetric media satisfying the relation 
e{uj)ii{uj) — 1 (of special interest in QCD), both with 
a sharp frequency cut-off on (^{lu) and with a behavior 
analogous to the one-absorption-frequency Sellmeir's for- 
mula for dielectrics. In Ref. for non- magnetic dis- 
persive media with spherical geometry, and based upon 
the Minkowski energy-momentum tensor, the Casimir 
surface force is worked out when the permittivity e(w) 
presents a step along the imaginary frequency axis. In 
the limit of perfect conductivity, a non-dispersive term is 
recovered which is in agreement with Boyer result |^,^ , 
while an attractive dispersive contribution is also found. 

In what follows we consider a simple model of a scalar 
field subject to frequency dependent local boundary con- 
ditions on the surface of a sphere. Thus our main goal 
is to establish a method for calculating the change of the 
Casimir energy of the field when the radius of the sphere 
is varied, in a situation where the boundary conditions 
impose a physical frequency cut-off Q. 

To this end we consider the very simple case of a scalar 
field whose modes corresponding to eigenfrequencies lu < 

are confined to the interior of a sphere of radius R, 
satisfying local homogeneous boundary conditions. 

On the other hand, we will assume that the boundary 
is completely transparent for those modes with uj > Q. 
Therefore, their contribution to the difference of Casimir 
energies for two different values of R will cancel out, 
no matter the regularization employed for its definition. 
Consequently, we will subtract these contributions, which 
amounts to a redefinition of the reference energy level in 
an i?-independent way. 

For the evaluation of the vacuum energy of the low fre- 
quency modes we will employ asymptotic expansions in 
an incomplete-C, summation technique, to be discussed in 
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the following. This approach will allow for the identifica- 
tion of the volume, surface, . . . , dominant terms in the 
Casimir energy. 

As a final exercise we will show that, with reasonable 
values for the frequency cut-off and the bubble radius, 
the amount of Casimir energy available in this model is 
comparable with the energy emitted during each cycle 
of a typical sonoluminesce experiment In view of 
the controversy existing on the subject, it would be of 
great interest to apply this method to a similar model 
for the case of the electromagnetic field in the presence 
of dielectric media, calculation which will be presented 
in a forthcoming paper p3]. 



II. THE MODEL AND ITS INCOMPLETE 
C-FUNCTION 

Let us consider a free scalar field in satisfying at the 
surface of a sphere of radius R, local boundary conditions 
which depend on the frequency lu of the field modes. 

We will make the assumption that the boundary is 
completely transparent for the modes of frequencies 
grater than a cut-off il, while for uj < il the modes satisfy 
Dirichlet boundary conditions, 



(A-f V^(f)-0, for r<i?. 



(1) 



being confined to the interior of the sphere. 

Writing '0w(r) = fi{r)Yf^{d, ip), we get for the radial 
function 



— + + — ]fi{r) = 0, for r < R, (2) 



dr^ r dr 



where the eigenfrequencies are determined by imposing 
the condition 



0. (3) 
The solutions regular at the origin are given by fi{r) = 



— Juiz), with V — I + 1/2 and z = w,y„r/c, where the 
2z 



eigenfrequencies are 



(4) 



being j^.„ the n-th zero of the Bessel function J^{z). 



i?-independent subtraction, which is nothing but a re- 
definition of the zero energy level. 

Therefore, we should evaluate the (finite) sum 



E{R)= ^"Y. o^'^-^"' 



(5) 



u=l/2 n=l 



where is the number of positive zeroes of Ji,(z) less 
than or equal to x = flR/c, the factor 2 j/ = 2 Z + 1 is the 
eigenvalue degeneracy, and vq is the maximum value of 
ly for which N^, > 1. 

We are interested in an analytic, rather than numeric, 
evaluation of eq. Although this is a finite sum, we 
will employ a summation method based on the evaluation 
of an incomplete (^-function, an approach which could be 
applied in more complex situations. We can employ the 
following representation: 



n=l 



(6) 



where the sum in the right hand side obviously exists for 
any s G C []. 

Since Juiz), for // > — 1, has only real zeros, and its 
non- vanishing zeros are all simple [ p4[ , we can employ 
the Cauchy theorem to represent the sum in the r.h.s. of 
(0) as an integral on the complex plane. 



n=l 



c JAz) 



(7) 



where the curve C encircles counterclockwise the N^, first 
positive zeros of Ju{z). 

For 5R(s) large enough, the contour C can be deformed 
into two straight vertical lines, one crossing the horizontal 
axis at 3fJ(z) = x and the other at 3fJ(z) = 0+. Indeed, 
expressing the integrand in terms of the modified Bessel 
function 



(8) 



(valid for — tt < arg(ii') < 7r/2), and taking into account 
its asymptotic behavior for large arguments [ pst , it is 
easily seen that, for < a; 7^ j^,„, Vn, the integral 



C,u{s,x) 



1 fX + lCO T' ( ~\ 

I ^-s'JjAz)^ 



2m 



Ju{z) 



(9) 



We will be interested in the difference between the vac- 
uum energies of two situations differing in the value of R. 
Then, we can disregard the contributions of those modes 
with frequencies ld > VL because, being independent of 
the position of the boundary, their contributions cancel 
out (whatever the regularization employed in defining the 
vacuum energy would be). This simply amounts to an 



^Notice that the sum in the r.h.s. of eq. evaluated at 
s = gives A*'^, the number of eigenfrequencies contributing 
to the Casimir energy of the field (after the subtraction made 
to define it) for a given value of the angular momentum / = 
V - 1/2. 
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converges absolutely and uniformly in the open half-line 
s > 1, from which it can be meromorphically extended 
to the whole complex s-plane. 
Therefore, for s > 1, 



E 



7-" 



C.(s,0+)-C.(s,x). 



(10) 



And, since the left hand side of (10) is holomorphic in s. 



the singularities of (^^{s^x) must be independent of x. 
On the other hand, for y > p^]. 



ly {~y -ix)^e ^'^"I^ (y + ix) 

Iv {y + ix) = (/^ {y - ix))* . 
So, for real s > 1 we can write 

as,x) = 



(11) 



3? ■ 



-^77(^ + 1) 



'-'-'ftlel^. 



(12) 



In order to construct the analytic extension of (^^ (5, x) 
to s ~ — 1, we subtract and add to the integrand in (|l^ ) 
the first few terms obtained from the uniform asymptotic 
(Debye) expansion |^ of the Bessel functions, 



I'A^t) 



-DAt)+Oiiy-^), 



Iu{vt) 
where 

D^{t) = iyD^^\t) + D^°\t) + i^-^D^-^\t) 



(13) 



(14) 



t 



2(1 + ^') 8^.(1 + ^2)1 



valid for large v with fixed t. We will see that this ap- 
proximation allows for the identification of the volume, 
surface, . . . contributions to the vacuum energy. 

Changing the integration variable in eq. (12h to t = 
z[y — i), with z = xjv > 0, we get 

C,v{s,x) = 



t-{s + l) d(ln/,(z.t)) 



dt 



(15) 



So, we must consider the integral 



— ^,.d(ln/.(.^))^^ 

— iz 



dt 



t-" D^{t) dt 



dt 



(16) 



The second integral in the right hand side of this equa- 
tion converges for s > — 2, since the integrand can be 
estimated by means of the next {0{i'^^)) term in the 
Debye expansion (eq. ([l3|)), which behaves as 0{t~^) for 
large \t\. It can be numerically evaluated at s = — 1. This 
will not be done in this paper. 

In the following, we will consider only the first inte- 
gral in the right hand side of (16), retaining the first few 
terms of its expansion in powers of consistent with 
the approximation made in eq. ([l^) . 

Notice that the integrand is an algebraic function, hav- 
ing singularities a.t t = 0,±i, and behaving as ©(t'^) for 
large \t\. This integral converges absolutely and uni- 
formly for s > 1, where it defines an analytic function 
to be meromophically extended to the region of interest 
of the parameter s. As we will see, this extension re- 
veals the singularities of C(s, x) as simple poles, whose 
residues are independent of x (a necessary condition to 
give a finite result in (10) for any s). 

In fact, by virtue of the analyticity of the integrand, 
we can change the path of integration to write 



t-' D^{t)dt 




dt 



t 

1 

+ 27^ 



it-t-^ 



2(l + t') 8^.(1 + ^2)^ 



1 

2t 



1 
2^2 



Svt^ 
1 

2t 



dt- 



8iyt^ 



dt. 



(17) 



The first integral in the r.h.s. of eq. (pT[), containing the 
whole dependence on x = vz, is holomorphic in s and 
can be directly evaluated at the required value of this 
parameter. On the half-line (l,oo) we have subtracted 
and added the first terms in the series expansion of Di, (t) 
for large t, thus making the second integral to converge 
for s > —2. The third one must be evaluated for s > 1 
and then analytically continued to the relevant values of 
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s. This can be exactly done, its contribution to (^i,(s,x) 
in eq. (|l^) being the real part of 



--TT + 

e ^ 

1/1+'^ 8 TT 



1 - s 



4z/ l-4i^2 

^ — ; 

s 1 + s 



(18) 



This expression has simple poles at s ~ 0,±1, which are 
the only singularities of (,j{s,x) for di{s) > —2. Notice 
that the residues of Ci^{s, x) are independent of x, 



ResC,,{s,x)\s=i = - , 

TT 



ResCi,{s,x)\s=o = 0, 

1-4^.2 



(19) 



ResCi^{s,x)\s=- 



8n 



and in agreement with the results in |^ (where Cv{s, 0"'") 
is studied). 

For example, for ^^{s,x) around s = — 1 and for ly < x 
(which will be needed in Section IV to evaluate the vac- 
uum energy), one straightforwardly obtains the Laurent 
expansion 



C„{s,x > v) 



Stt (1 + s) 
log (^0 + log (^z + \J -\ 



V 

47r 



-zv/.^-l 



24 TT (z2- 1)3/2 37, 

+0(.s + l). 



(20) 



with fixed z = xjv^X. 

On the other hand, for a; ^ 0+ a similar calculation 
leads to 



C.(5,a; = 0+) 



1-4;/" 
Stt (1 + s) 

1 



47r 



Stt 



^^logn 
Stt ^V2/ 

+ 0(s + l). 



(21) 



In the following Section we will evaluate, as a function 
of i/, the number of modes contributing in eq. (|l^, and 
in Section IV , their contributions to the vacuum energy. 



III. THE NUMBER OF CONTRIBUTING MODES 

In this Section we address ourselves to the determi- 
nation of vq in (P), the maximum value of v for which 



Njj > 1. Although in the simple case under study the 
zeros of Jv{w) are well known |2^, we prefer to estab- 
lish a criterium which can be applied in more general 
situations. 

First, notice that 



n=l 



= [Us,0+)~Us,x)]\^^, (22) 



s=0 



is a discontinuous function of x, having a step of height 
1 at each positive zero of the Bessel function J^, (w) . 
Then, can be determined from the condition 



N,^,{x)=N,,ij,„^i+0)^l, 



(23) 



with A^^„(>o,i -0) = 0. 

Taking into account eq. (^9|) and the fact that the sec- 
ond and third integrals in the r.h.s. of eq. ( p^ are real, 
it is straightforward to obtain from eqs. (p"^^p^ that 

C.(s = 0,a;) = -^-i + 
+^{"-^ (^yTTF*^-log(l + y/l + e-''' z^)^ 



4 TT 



z(2-|-3z2) 



24t/7r(l + e-»'^z2)3/2 



(24) 



where we have taken z = x/i^ 1. In particular, for 
X 0+, 

C,(s = 0,x = 0+) = -^-i+O(i.-2), (25) 

in coincidence with [p6| . 

Now, taking the difference between ( p5| ) and ( p^ we 
get a smooth approximation, N^(x) -1-0(1^"^), to the step 
function Ny{x) in (p2[). It is easily seen that, for v > x, 
Ny{x) = while, for < x, we have 

N^{x) = - (Vz^ - 1 - arctan(v/z2 - 1)) - 



2-F3z2 



(26) 



4 24i/7r(z2 - l)='/2' 

with z = xjv. 

Let us now determine the value vq for which 

iV.„(x) = l/2. 

To this end, we propose an expansion of the form 

r~2 7 -1/3 , -3/3 , -5/3\ 



(27) 



(28) 
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which makes sense for vq ^ 1 and zq = x/pq ~ 1- Re- 
placing in (|26| ) and imposing (p7|), the coefficients Sk can 

— 1/3 

be determined order by order in Vq , to get 



x = vo + 1.857 i^o^^ + 1.034 i/g + ©(i^, 



) 



or, inverting this development, 

vq = x-1M7x^^^ + 0.1155 x-^^^ + 0{x-^). 



(29) 



(30) 



(Notice that vq < x.) 

Equation (^) is in excellent agreement with the ex- 
pression of the ffi'st non- vanishing zero of Ji,g(w), for 
a large order vo ||2^]: j^^^i — 



1.8557 1', 



1.03315;/, 



-1/3 



1/3 




OK"')- 



IV. THE DOMINANT CONTRIBUTIONS TO 
THE VACUUM ENERGY 

In this Section we evaluate the first contributions to 
the vacuum energy obtained from the Debye expansion 
employed in Section ^ As we will see, this allows to 
determine the volume, surface, . . . , terms in the Casimir 
energy of the scalar field. 

According to the results in the previous Section, we 
are interested in the Laurent expansion of C^(s,a; > v) 
and Ci/(s, X = 0+) around s = —1, given in eqs. ( pO| ) and 
( ^ ) . As already remarked, the singular parts cancel out 
in the difference on the r.h.s. of eq. ( p^ (see eq. (p^). 
For the difference of the finite parts we get 

zF{x,v)= [C,(s,G+)-C.(s,x)]^_^ = 



2 

^ (z V-l + - log(z + V-1 + z^)) 

2 TT V / 



3z- 8z^ 



247r(-l + z2)2 



-Oi 



l0g(z + V-l + z2) 



(31) 



This is a good approximation as long as ^ 3> 1 and z 
x/v > 1. 

Our aim is now to evaluate the sum in eq. (|^), 

E{R ^xc/VL) ^ 



^ y: .[c.(-i,-)-c.(-i,o+)] 



(32) 



The function F{x,i/) is non-negative and has a pro- 
nounced maximum a,t v ~ x/2 (i.e. z « 2). Thus, the 
use of the approximation in eq. (31 ) is consistent if x 3> 1. 

From eq. (31), it is not difficult to sec that the succes- 
sive terms in the Euler - Maclaurin summation formula 
11, 

V F{x,iy) = / F{x,v)dv+ 



v=l/2 



+ -{F{x,v^) + F{xM2)} + 



(33) 



12 



9„ [F{x, y)] 



+ 



v=l/2 



are of increasing order in x^^ . So, retaining the first few 
terms consistent with the approximation made, we get 
for the vacuum energy 

E{R) ^ 

nn 

-(V-l+^g(-5+2^g)) 2^-3 log(^o+V^mi) 
247rz^ ' 24 tt z% 



^2 / -y-l + ^g , 7r+3 log(zo + V-l+zg) 
X I I TtTT3 



-(-27+17 2g) _^ 6 7r-10zg+69 logCzp + ^-l+zg) 
48 20 V-l+-^o 



48 7rzg 



-O(xO), 



(34) 



where zq = xjvi^^ x — Ril/c and vq is given in eq. (|3^). 
In this equation one can recognize volume, surface and 
curvature contributions to E{R). Notice that we could 
have retained any number of terms of the asymptotic 
expansion in eq. (p^), and performed the same steps as 
in the present calculation, to get the Casimir energy to 
any order in x^^ . 

Finally replacing in eq. (|^ zq x/vn and z/q by its 
expression we get 



E{R) ^hn 



1277 12 



0.1343x3 +0{x) 



(35) 



where one can see that the volume and surface terms 
are dominant for a; >> 1. Notice that also non-integer 
powers of the radius R appear as a consequence of the 
relation between i^o and x, eq. (pO|). 



iy=l/2 



with vq given in (| 



V. CONCLUSIONS 



In equations ( |34|]35| ) we have derived the dominant con- 
tributions to the vacuum energy of a scalar field in a 
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model with a frequency dependent boundary condition, 
consisting in the confinement of the modes with low fre- 
quency (up to a physical cut-off fl) to the interior of a 
sphere of radius R. 

These modes are subject to Dirichlet boundary condi- 
tions at the surface of the sphere, while those with fre- 
quency higher than Q are free, being the boundary com- 
pletely transparent to them. This characteristic of the 
model allows for the subtraction of the contribution of 
the high frequency modes to the vacuum energy, which 
amounts (independently of the regularization employed 
to define it) to an i?-independent redefinition of the zero 
energy level, having no consequences on the evaluation 
of energy differences. 

In so doing, we have represented the sum over the 
eigenfrequencies up to the cut-off Q in terms of an in- 
complete ^-function associated with the Laplacian oper- 
ator in the sphere with Dirichlet boundary conditions 
(see eq. (|l0|)). The function Ci^j^) i^i 6q. (^2|) is well 
defined only for ^{s) > 1. So, it was analytically contin- 
ued from s > 1 to the relevant values of this parameter 
(s = 0, needed to evaluate the maximum angular mo- 
mentum Iq = vq — 1/2 giving rise to eigenfrequencies less 
than or equal to Q, and s = —1, necessary to evalu- 
ate the contribution to the Casimir energy of the modes 
with angular momentimr I = v — 1/2) by approximat- 
ing the behavior of the integrand in eq. (|l2) employing 
the Debye asymptotic expansion of the modified Bessel 
functions appearing in its expression. 

This procedure has lead to a meromorphic function 
having simple poles with (exactly evaluated) cut-off in- 
dependent residues (see eq. (pj|)), a necessary condition 
to have a finite result for the sum in eq. (10). 



The finite part of Ci^^^) has been evaluated up to 
terms of a given order in u~^. Although we have re- 
tained only the first terms in this asymptotic expansion, 
one can follow the same steps to determine the finite part 
of <^(s, x) at s = —1 or up to any required precision in 

Finally, the application of the Euler - Maclaurin sum- 
mation formula has lead to an expression for the Casimir 
energy of the model in which one can recognize volume, 
surface and curvature contributions (see eq. (^J)). 

For a cut-off corresponding to x = RVl/c » 1, the 
dominant terms in the vacuum energy, eq. (^Sf), are 
proportional to the volume {V = 47ri?'^/3) and area 
(S" = 47r_R^) of the sphere. 



V ■ 



16 TT^ 



12 7r2 c2 



(36) 



with ^ = -7r/4. 

It is worthwile to remark that, for a similar model 
where the low frequency modes of the scalar field are 
subject to Neumann (rather than Dirichlet) boundary 
condition, we get the same expression for the dominant 
terms with ^ = -|-7r/4. 

These two dominant terms are in complete agreement 
with those obtained from the expansion of the density of 



states in powers of the inverse wavelength^ 00. The 
relation between i/q and x, eq. (|30| ) (or, equivalently, the 
expression of the first zero of the Bessel function J^^ (w) 
in terms of the order vq ) introduces also non-integer 
powers of the radius R (see eq. (|35|)). 

As a final exercise, we can use eq. ( ^5| ) in a very 
schematic model pretending to mimic the phenomenon of 
sonoluminescence. We will adopt the values of the radius 
and emitted energy corresponding to a typical sonolumi- 
nescent bubble [Q, and estimate the cut-off Q needed to 
produce this amount of energy. To this end, we will sim- 
ply take the difference of the low frequency contribution 
to the Casimir energy of the scalar field for two different 
values of the bubble radius. 

If the bubble collapses from an initial radius R — 4 x 
10^^ m to a final radius of one tenth this value, and the 
emitted energy is £ — 1.2 x 10"^^ Joule, by imposing the 
equality 



— {E{R) - E{R/10)) = 1.516 X 10^ 
Tic 



(38) 



490, 



and taking into account eq. (|3g) , it follows that x 
justifying the use of the approximation obtained. 

This implies that ft = 3.675 x 10^^ 1/sec, which cor- 
responds to a cut-off in wavelengths in the ultraviolet of 



A = 5.129 X 10" 



5129 A, not far from the region 



where the refraction index of water becomes essentially 1 
psf . This strongly suggests to consider a similar model 
for the case of the electromagnetic field in the presence 
of dielectric media, calculation which will be presented 
elsewhere |^ . 
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^Indeed, for scalar fields subject to local homogeneous 
boundary conditions the density of states is modified by fi- 
nite volume efi'ects . The first correction in the asymptotic 
expansion for large wavelength k is given by 



E 



V 



(27r)3 



+ s i 



(2 7r)3fc 



+ . 



(37) 



where the coefficient ^ takes the value ^ = — 7'"/4, ^ = +7r/4 
and 5 = -|-7r/4 for Dirichlet, Neumann and Robin boundary 
conditions respectively. Then, introducing the dispersion re- 
lation (jj{k) = ck and a cut-off in the wavelength given by 
K = Q./c, it is easy to get eq. (^) for the vacuum energy. 

For junction boundary conditions, as those appearing in 
problems with dielectrics, the coefficient ^ has more involved 
expressions, which can be difficult to evaluate 11 1|. 
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